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Helical magnetic fields are injected into the cosmic medium during cosmological baryogenesis and 
can potentially provide a useful probe of the early universe. We construct a model to study the 
injection process during a first order phase transition and to determine the power spectra of the 
injected magnetic field. By Monte Carlo simulations we evaluate the Fourier space symmetric and 
helical power spectra of the magnetic field at the time the phase transition completes. The spectra 
are peaked at the scale given by the inverse size of bubbles at percolation and with a comparable 
width. These injected magnetic fields set the initial conditions for further cosmological magneto- 
hydrodynamical evolution. 



Our study of the universe relies on relics left-over from 
early cosmological epochs. The cosmic microwave back- 
ground brings information from the epoch when atoms 
formed, the light elemental abundances from the epoch 
when nuclei formed. Similarly, the electroweak phase 
transition may mark the epoch when a net baryon num- 
ber was generated, or when net lepton number was con- 
verted into net baryon number, and this coincides with 
the generation of magnetic fields. The present paper is 
based on the hypothesis that primordial magnetic fields 
will inform us of the epoch when a net amount of baryons 
first formed, the so-called epoch of "baryogenesis" [H, 0] . 

The question of whether a primordial magnetic field ex- 
ists is often raised in connection with the magnetic fields 
observed in galaxies and clusters of galaxies with strength 
~ fiG and kpc-Mpc coherence scale. There is consider- 
able debate whether primordial fields arc essential to the 
generation of galactic fields, and what properties of the 
primordial field are necessary to turn them into observed 
magnetic structures. The arguments involve the coher- 
ence and amplitude of observed magnetic fields, the effi- 
ciency of galactic dynamos, the turnover time scales asso- 
ciated with galactic dynamics, especially with the earli- 
est known galaxies containing magnetic structures, astro- 
physical sources e.g. active galactic nuclei that may spew 
out magnetic fields, and the generation of large scale seed 
fields by the Biermann battery. We shall by-pass these 
issues since, in our view, a primordial magnetic field is of 
interest in itself, whether or not it is responsible for the 
observed magnetic fields in galaxies. If there are strongly 
motivated early universe scenarios, based on reasonably 
well-established particle physics, that lead to the genera- 
tion of magnetic fields, they provide good reason to study 
and to look for these structures in cosmological data. 

The connection between baryon number production 
and primordial magnetic fields can be understood intu- 
itively in the following way. Baryon number violation 
in the standard model of the electroweak interactions is 
made possible due to a quantum anomaly. As a physical 
process, baryon number is generated when many differ- 
ent particles come together to form an object called a 



"sphaleron" [U, 3 , which then decays into a final state 
with baryon number that is different from that of the 
initial state. The sphaleron, and its deformations, are 
made up of scalar and electroweak gauge fields and may 
also be viewed as being an unstable bound state of an 
electroweak magnetic monopole and an antimonopolc, 
with an electroweak Z— string confining them @, 0j- The 
monopole and antimonopole attract each other by the 
Coulomb force and the confining Z— string also pulls 
them together, and a static solution is in general not 
possible. In the sphaleron though, the antimonopolc has 
a twist in field space relative to the monopole which pre- 
vents the monopole and antimonopolc from annihilating, 
and allows for the existence of a static solution 0, lf|. 
The presence of magnetic charges also explains the large 
magnetic moment of the sphaleron calculated in Ref. [J] . 
Twisted solutions similar to the sphaleron also occur in 
the context of kinks in one spatial dimension and can 
lead to novel static phases containing a lattice of kinks 
and antikinks [![. 

Once we appreciate that baryon number violation pro- 
cesses have intermediate states that consist of monopolc- 
antimonopole pairs, it is not hard to see that magnetic 
fields must be produced when a sphaleron decays. But 
there is a further "twist" to this connection. The twist in 
the fields that prevents the monopole and antimonopole 
from annihilating is not stabilized, and the monopole can 
untwist and annihilate the antimonopole. The instability 
of the system causes the sphaleron to decay and radiate 
away its energy, releasing magnetic fields in the process 
@. The decay has been studied numerically and a very 
interesting feature emerges at late times. The released 
magnetic fields inherit the twist of the sphalerons and is 
measured by the magnetic hclicity integral 



n = / d 3 xA ■ B 



(1) 



At late times, the magnetic field evolves such that the 
magnetic hclicity is conserved. The conservation of 
magnetic hclicity is familiar in magneto hydro-dynamics 
(MHD) in plasmas with high electrical conductivity. Yet 
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in Q there is no external plasma; the only charges in 
the system are those resulting from the decay of the 
sphaleron itself. Remarkably, magnetic helicity is still 
conserved during sphaleron decay. 

To summarize, baryon number violating processes in 
the electroweak model occur via the production and sub- 
sequent decay of sphalcrons. Each sphaleron produces 
helical magnetic fields when it decays and the helicity of 
the magnetic field is conserved. Now, since the produc- 
tion of each baryon gives a certain amount of magnetic 
helicity, the cosmological baryon number density, , can 
be related to the magnetic helicity density (h) 



h w — rib 



(2) 



where the minus sign requires more detailed considera- 
tions and indicates that the primordial magnetic field 
is left-handed. (In writing Eq. ([2]) we are assuming that 
the cross-helicity between the magnetic field produced 
by different sphalcrons averages out to zero.) In princi- 
ple, there could be a numerical pre-factor on the right- 
hand side and, in fact, early estimates suggested it to be 
~ a -1 where a — 1/137 is the fine structure constant. 
However, explicit numerical evolution of the sphaleron 
suggests that @ may be a better estimate [{|. Numer- 
ically, the cosmic baryon number density is ~ 10 /cm 3 
at the electroweak epoch. 

The rough equality of magnetic helicity and baryon 
number gives us an estimate for the integral in Eq. ([1]) 
but does not provide much information about the char- 
acteristics of the magnetic field itself. Of interest are the 
two point correlation functions of the magnetic field. As- 
suming statistically homogeneous and isotropic magnetic 
fields, the spatial correlator at a fixed time can be written 
as 

CyW^Cx^x + r)) 

= - J d 3 x Bi(x)B J (x + r) 

= M N {r)Pij + M L {r)fifj + e ljk f k M H (r) (3) 

where, V is the integration volume, k = 1,2, 3, r = 
|r|, f = r/r, e,ij k is the Levi-Civita tensor, P^ is the 
traceless projection tensor, 



Pi 



(4) 



and Mjv, Ml and Mh denote the "normal", "longitu- 
dinal" and "helical" correlation functions respectively. 
Also note that the helical term in Eq. has a factor 
of f k and not r k in it, as is sometimes used. The cor- 
relation functions will also depend on time but here we 
are only interested in the injected field at the end of the 
electroweak phase transition. So it is to be understood 
that the correlation functions are all evaluated at t = t ew 
where t ew denotes the time at which the phase transition 
is complete. 

The normal and longitudinal correlation functions are 
not independent since the magnetic field is divergence- 
less, V ■ B = 0, and the relation between Mm and Ml 



d -^ = 2 -(M N -M L ) 
dr r 



(5) 



It is more conventional to work in Fourier space, where 
the correlator can be written as 

<6*(k)&,(k')) = 

[S(k) Pij + ^e ^ ■ ^ fc ^ A(fc)] (2^) 3 <5( 3 >(k - k') (6) 

with no longitudinal term (proportional to kikj) since 
the magnetic field is divergenceless. We have denoted the 
Fourier components of B(x) by b(k) with the convention 



6i(k) - J d 3 xBi(x)e 



?kx 



and the momentum space projection tensor is 



Pij 



Si 



(7) 



(8) 



By explicitly Fourier transforming dj (r) it is possible 
to derive 



S(k) = Ait I dr 



Mr 



sin(fcr) cos(fcr) 



(At) 



M, 



N 



sin(fcr) sin(fcr) cos(fcr) 



{krf 



{kr)< 



(9) 



Further simplification occurs on using the divergenceless 
condition in Eq. (0, 



cos(At) )Ml (10) 



where, to omit the boundary term, we have assumed 
r 2 Mi(r) — > as r — > oo. While Eq. (fTOj) is simpler, we 
still use Eq. ((HI) in our numerical work because r 2 M^(r) 
is not negligible at the scale of our simulation box. Sim- 
ilarly we have the relation 

4tt f°° , /sin(Ar) „ A 
A(k) = — J dr r f ± ' - cos(kr)j M H (11) 

The goal of this paper is to develop a model for the 
generation of magnetic fields during baryogencsis (Sec. [J), 
to Monte Carlo the model (Sec. [TTJ) and to obtain the 
correlation functions of the magnetic field fSec. HII[) . Our 
final results will be the functions Mj\r(r), Ml(t), Mn(r), 
S(k) and A(k). The spectral properties of the injected 
field should be useful to study their evolution. Wc discuss 
limitations of our model and future prospects in Sec. IIVI 



I. MODEL 

The relation between baryogenesis and magnetic fields 
exists as long as there is anomolous baryon number vi- 
olation. Even in leptogenesis scenarios, sphaleron pro- 
cesses arc required to convert lepton number into baryon 
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number [10( , and these conversions will produce magnetic 
fields. The power spectra of fields produced during lepto- 
genesis, however, may be different from those produced 
during baryogenesis as considered here because, unlike 
baryogenesis, leptogenesis does not rely on a first order 
phase transition for departures from thermal equilibrium. 




We assume that baryogenesis occurred at the elec- 
troweak scale and the necessary departures from ther- 
mal equilibrium are provided by a phase transition that 
was strong enough. The electrowcak symmetry is broken 
within bubbles which then grow, merge and eventually fill 
all of space, thus completing the phase transition. During 
this process, baryon number changes can occur relatively 
freely outside the bubbles since there is little energy cost 
associated with the sphaleron in the symmetric phase, 
but baryon number changes are highly suppressed within 
the bubbles since the sphaleron has mass ~ mw /ot in the 
symmetry broken phase. (Here mw ~ 100 GeV is the 
mass of the electroweak W gauge boson.) Also, in the 
symmetric phase that is outside the bubbles, the electro- 
magnetic magnetic field is just a linear combination of 
the electroweak magnetic fields, and has no special sig- 
nificance. For example, all the electoweak gauge fields 
are massless, and we expect rapid interactions to main- 
tain equilibrium between the different degrees of free- 
dom. Inside the bubbles, the electromagnetic field is the 
only massless gauge field and its evolution is described by 
Maxwell's equations. However, there are no sphaleron 
transitions inside the bubble and so no magnetic fields 
are generated there. It is only the sphaleron transitions 
occuring in a thin region right around the bubble wall 
which produce the magnetic fields that are then captured 
by the growing bubble. It is this magnetic field, gener- 
ated by sphaleron transitions at the surface of bubbles, 
that is of interest to us. 



A schematic picture of our model is shown in Fig. [T] 
Essentially, there arc bubbles of the broken symmetry 
phase that nucleate and grow, and sphalcrons are explo- 
sions that occur on the surfaces of these bubbles and blow 
out magnetic fields into the environment. We would like 
to find the correlation functions of the magnetic fields 
left-over after the phase transition has completed. Our 
picture of the phase transition is similar to that of cos- 
mological large-scale structure formation in which cosmic 
voids nucleate and grow. Supernovae or other astrophys- 
ical activity occurs on the surfaces and intersections of 
the voids and expels magnetic fields and other elements 
into the cosmic environment. 
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FIG. 1. Schematic representation of the electroweak first or- 
der phase transition via growing bubbles (circles), and of the 
sphaleron explosion-like events that occur on the surfaces of 
the bubbles. The discs represent the magnetic fields from the 
explosion event, with earlier explosions having had more time 
to grow out further. 



a sphaleron explosion as 



B r = 



a cos 9 



(a 2 + r 2 )3/2 



2)5/2 



(12) 



-e- r / a sin< 



where (r, 9, 0) are spherical coordinates centered at the 
location of the sphaleron and the z-axis is chosen so that 
there is azimuthal symmetry. (Plots of two sections of 
the magnetic field are shown in Fig. [2]) The r, 9 com- 
ponents of the magnetic field are chosen to be approxi- 
mately those of a current-carrying circular loop of wire 
of radius a, and with current proportional to 1/a. (The 
approximation in Eqs. (|12p for the field of a circular loop 
of wire is taken from Sec. 5.5 of Ref. [H|.) The azimuthal 
component of the magnetic field is chosen so that it is lo- 
calized in a region of size ~ a. The size of the whole sys- 
tem is chosen to grow with time to model the exploding 
sphaleron. Assuming that the magnetic fields generated 
by the sphaleron expand out at the speed of light, we 
can take a(t) = t — to where to is the epoch at which the 
sphaleron transition occurred. 

The scalings with a are important since the magnetic 
helicity in the aftermath of a sphaleron decay stays con- 
stant, as we now show explicitly. We first construct the 
gauge potential for the azimuthal component of the mag- 
netic field 
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(13) 



The only element that is missing from the picture so far 
is a model for the "explosion" that expels the magnetic 
field. In Ref. the decay of the sphaleron was studied 
numerically, with the result that magnetic helicity stayed 
constant at late times, and the energy density spread out 
as time progressed. We shall model the magnetic field in 



In other words, 



Now we define 



V x (4 azim) 0) 



b = bw+b< 2 > 



(14) 
(15) 
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FIG. 2. Plots of the magnetic field vectors in the xz-plane 
(top) and xy-plane (bottom) with a = 1 in Eq. (|12p . 



where, 



B (1) = B r f + B e 



B< 2 > = B^4> 

and, correspondingly, 

A= A«+A^ 



(16) 



(17) 



(18) 



so that BW = V x AW. Then since the helicities in B^ 
and B' 2 ) are individually zero, the total helicity is given 
by the cross terms 

H = J d 3 xA^ ■ B( 2 ) + J d 3 xA^ ■ B« (19) 

An integration by parts is now used to obtain 

H = 2 I d 3 xA ( - 2) ■ B (1) (20) 



The boundary term vanishes because the fields are local- 
ized. 

Next we insert the expressions for the gauge potential 
and magnetic field as given in Eqs. (fT"2"]) , ([13]) , A simple 
change of variables, u = r/a, shows that the magnetic 
helicity of one source is 



H= —0.57 = 4.: 



(21) 



for any value of a. Therefore the magnetic helicity is 
conserved for our choice of magnetic fields and matches 
the conservation seen in sphaleron decay 0] . 

The choice of factors of a in Eq. (Ti~2"]) also ensures that 
the relative energy in each of the three components of 
the magnetic field stays fixed, while the net energy in 
the magnetic field decays as 1/a. The decay of the mag- 
netic energy is a necessary consequence of the conserva- 
tion of helicity because the energy density is B 2 while the 
helicity density is A • B and hence, simply by counting 
dimensions, the total energy is the total helicity divided 
by a length scale. Since the only length scale in the prob- 
lem is a and the total helicity remains constant, the total 
energy must decay as 1/a. 



The magnetic field in Eq. ([12]) is axially symmetric 
about the chosen z-axis. However, different sphalcrons 
will produce magnetic fields that are azimuthally sym- 
metric with respect to different axes. So a sphaleron 
is described by its location as well as its orientation. 
We will assume that the orientations of the sphalerons 
are isotropically distributed and, for a given sphaleron, 
choose the orientation from a uniform distribution on 
the two-sphere i.e. in spherical coordinates, cos# and <p 
are chosen from a uniform distribution over the intervals 
(—1,1) and (0, 27r) respectively. In principle, the inter- 
action of electroweak fields with the bubble wall could 
result in a preferential orientation of the sphaleron (e.g. 
normal to the bubble wall) but we shall disregard this 
possibility in the present paper. 

We do not claim that an electroweak sphaleron pro- 
duces the magnetic field in Eq. (fT2"]) when it decays. In- 
stead, (|12p is a convenient choice for the magnetic field 
and has the following desirable properties: (i) The field 
is smooth and divergenceless. (ii) The magnetic helic- 
ity is independent of a. (iii) The relative energy in all 
three components of the magnetic field is independent of 
a. Furthermore, we arc only interested in the correlation 
functions for the magnetic field at large separations and 
hope that these are not sensitive to the exact form of the 
model we choose for the sphaleron's magnetic field. To 
put this work on a firmer footing, it will be necessary to 
study the process of sphaleron decay more carefully and 
to devise a more accurate model for the magnetic fields 
produced. 
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II. MONTE CARLO SIMULATION 



The procedure we follow is to throw bubble sites ran- 
domly with uniform distribution within our simulation 
volume. Successful bubble sites are those that lie outside 
of all existing bubbles. Then the bubbles grow at speed 
Vb. The bubble growth velocity depends on the ambient 
plasma and can be much smaller than the speed of light 
(e.g. pjjl, [H| , [3). However, to keep the number of 
parameters to a minimum in our simulations, we took 
Vb = c = 1. As the bubbles grow, we randomly nucleate 
sphalcrons on the bubble surfaces with angular density 
1/A S where A s denotes the average area per sphalcron. 
Here also we take care to eliminate sphalcron sites that lie 
within pre-existing bubbles. Each sphaleron is described 
by its time of nucleation, location, as well as its orienta- 
tion. The sphaleron nucleation time enters the factor a 
in Eq. (fl"2|). while the spatial location sets the origin for 
the axes, and the orientation of the sphalcron fixes the 
direction of the local z-axis. With time, the radius of the 
magnetic field grows with velocity v m which, for conve- 
nience, is also taken to be the speed of light, v m = c = 1. 
After the phase transition is complete - no more bub- 
bles nucleate since almost all the simulation volume is 
occupied by pre-existing bubbles - we find the magnetic 
field on a lattice within our simulation volume due to all 
sphalcrons in our simulation box. Wc then compute the 
correlation function, CV,(rf) (Eq. which has 3x3x3 
components because of the 2 free indices i,j and the 3 
choices for the direction of f. The simulation is run many 
times with different seeds for the random number gen- 
erator and ensemble averages are calculated. We have 
explicitly checked that the correlators are of the form in 
Eq. ([3]). Then the spectral functions are found as linear 
combinations of the CV,(rf), 



f 

M N (r) = ^P ij C ij {rf) 

f 

M H (r) = -^eijkhCijirf) 



(22) 
(23) 
(24) 



where the sum is over the 3 directions: f € (x, y, z). 

The simulation takes care that bubbles can only nu- 
cleate in the false vacuum region, that is, outside every 
other bubble. In practice, a certain number, of bub- 
bles are thrown down at every time step, and any bubbles 
that lie within pre-existing bubbles are rejected. Next we 
want to locate sphalcron events on the surfaces of existing 
bubbles. We let A s denote the mean area occupied by a 
sphaleron on the surface of the bubble. On the surface of 
a bubble of radius R, we throw 4irR 2 /A s sphaleron sites 
where we take A s = 32 so that there are a large number of 
sphalerons (order 10 5 ) in the simulation but still within 
computational limits. We reject those sites that lie within 
any other bubble. For large bubbles, the mean distance 




FIG. 3. The number of bubbles nucleated at time t for 3 
different nucleation rates. The nucleation rates correspond to 
throwing down 30 (black solid curve), 40 (dashed red) and 50 
(dotted blue) bubble sites per time step. The bubbles rapidly 
fill the simulation box and even with the lowest nucleation 
rate the phase transition is essentially complete by t = 25 i. e. 
25 time steps. 



between neighboring sphalcrons on the same bubble is 
« 2^/A s /{4ir) « 3. 

The bubble nucleation rate is chosen over a range that 
gives ~ 10 2 bubbles in the simulation. The sphaleron 
nucleation rate is kept fixed, while the bubble nucleation 
rate is taken to be 30, 40 and 50 bubbles per time step. 
The number of bubbles nucleated per time step for the 
three different nucleation rates are shown in Fig. |3] The 
error bars denote la fluctuations about the mean taken 
over 20 runs. In Fig. [4] we plot the number of sphalerons 
nucleated per time step, also including la error bars. 
This plot is equivalent to plotting the surface area sepa- 
rating the true vacuum and the false vacuum as a func- 
tion of time. 

A subtle point about the simulation is that we have 
nucleated bubbles within a box but then the bubbles 
subsequently grow beyond the box. The sphalerons are, 
however, nucleated only within the box at a fixed rate, 
rejecting those sphalerons that lie on the parts of the 
bubbles that are outside the box. So the magnetic field 
close to the boundaries of the box suffer from boundary 
artifacts due to the lack of sphalerons outside the box. 
Hence it is important that the magnetic field only be cal- 
culated in a sub-box that is smaller than the original box, 
at least by a margin that is larger than the size of the 
typical magnetic structure, given by a(t). In our sim- 
ulations, the box size was 144 lattice spacings and the 
sub-box size was 108 spacings, i.e. we excluded a bound- 
ary layer of 18 lattice spacings all around the box. As 
can be seen from Fig. 2J most sphalerons nucleated at 
t ~ 15 and had a ~ 10 when wc stopped the simulation 
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III. RESULTS 




FIG. 4. The number of sphalerons nucleated at a time t for 
the three different bubble nucleation rates as in Fig. [3] 



(t = 25). So a for most bubbles is less than the thickness 
of the excluded boundary layer. The other relevant sim- 
ulation parameters are the lattice spacing, dx = 1, and 
the time step, dt = 1. 

After the phase transition is complete, we know where 
all the sphalerons are located, their orientations, and 
also their sizes because we know the times at which the 
sphalerons exploded. We take the expansion speed of 
the magnetic fields produced by a sphaleron to be the 
speed of light. This determines the size, a(t), occurring 
in Eq. (|12p for every sphaleron. Then at every point on 
a sub-lattice we sum over the magnetic field due to every 
sphaleron. (This is the computationally expensive part 
of the code since it involves roughly 10 5 x (108) 3 ~ 10 11 
computations.) Once we know the magnetic field at each 
lattice site, we calculate spatial correlations by doing the 
volume integral in Eq. ([3]) and averaging over 20 en- 
sembles. Projections of the correlation functions as in 
Eq. (|24j) immediately give the normal, longitudinal and 
helical power spectra. The integrals in Eq. ©, (TTT1) fi- 
nally lead to the Fourier space power spectra, S(k) and 
A(k). 

Before proceeding to the numerical details and re- 
sults, we summarize the electroweak and cosmological 
parameters. The Hubble distance at the electroweak 
epoch is « 10 cm, while the thermal length scale 



is TL 



10 



-16 



cm. The sphalerons are exploding on 



the inverse electroweak mass length scale which is com- 
parable to the thermal length scale. The ejected mag- 
netic fields produced can spread out freely until the MHD 
frozen-in length scale at the electroweak epoch, Zf r0 zon ~ 
^/i cw /47rcr c « 10~ 8 cm, where a c ~ T cw /e 2 is the elec- 
trical conductivity of the plasma. The present baryon 
number density is n& ~ 10~ 7 cm -3 and at the electroweak 
epoch this corresponds to iib.cw ~ 10 37 cm -3 . 



In Fig. [J] we have already shown the nucleation of 
sphalerons as a function of time for each of the different 
bubble nucleation rates. (We do not vary the sphaleron 
nucleation rate.) The number of sphalerons nucleated at 
any time is proportional to the net surface area of the 
bubbles. As expected, the plot shows that the surface 
area, and hence the sphaleron nucleation rate, grows to a 
maximum and then decreases. If we increase the bubble 
nucleation rate, the rate of growth of the surface area 
is larger initially, but then the turning point is at earlier 
times because it is determined by the merging of bubbles. 

The plot of the sphaleron rate versus time in Fig.2]also 
tells us the size distribution of the expanding magnetic 
field distribution. In particular, lower bubble nucleation 
rates, as are relevant in strongly first order phase tran- 
sitions, lead to a sphaleron rate that is larger at later 
times. (The peak in Fig. 2] is shifted to the right.) So 
at some fixed late time, the sphaleron explosions have 
had less time to grow and there are a larger number of 
smaller sphaleron explosion remnants. The correlation 
functions should therefore be larger at small distances 
when the bubble nucleation rate is smaller. This fea- 
ture can be seen in Figs. [5] and [5] where we show the 
spatial correlation functions, Ml and M/v, versus r. In 
Fig. [7] we show the helical correlation function, Mjj (r) . 
The peak shifts to the left for smaller bubble nucleation 
rates in agreement with our observation above that more 
sphaleron explosions occur later if the bubble nucleation 
rate is small and have less time to grow. 

The fluctuations in the correlation functions, denoted 
by the error bars, are quite large. To further reduce them 
would require increasing the sphaleron rate and would 
increase the computational time. At present, each Monte 
Carlo with 20 runs takes about 10 days to run. 

The Fourier space correlation functions can be found 
using Eqs. fl§J, (fTTj) and are shown in Figs. [8] and |U The 
spectra are dominated by peaks at k ~ 0.05. This corre- 
sponds to a length scale I ~ k~ 1 ~ 20. From Fig. @] we 
see that most sphalerons were nucleated at t ~ 15 and 
these would primarily be on bubbles that are also of size 
~ 15, since most bubbles are nucleated at the earliest 
times. Hence the peak of the correlation is given by the 
sizes of the bubbles when bubbles started to percolate. 
The width of the peaks in Figs. [8] and U are Afc ~ 0.05 
and also given by the bubble sizes at percolation. Note 
that it would not be suitable to characterize the injected 
spectra by power law fits. 



IV. CONCLUSIONS AND FUTURE 
DIRECTIONS 

Following the general scenarios discussed in Refs. [l], Q 
we have proposed a concrete model for the generation 
of helical magnetic fields during baryogenesis at a phase 
transition. The model takes into account magnetic field 
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FIG. 5. Ml(j-) for different rate of bubble nucleation. Solid 
(black), dashed (red), and dotted (blue) curves correspond to 
30, 40 and 50 bubbles nucleated at every time step. Fluctu- 
ations are only shown for the run with 30 bubbles nucleated 
per time step. 




FIG. 6. Mjv(r). Plots are made following the scheme of Fig. [5] 

generation due to baryon number violating processes oc- 
curring on bubble walls. By Monte Carlo simulations, we 
have evaluated correlation functions of the injected mag- 
netic field on completion of the phase transition. The 
Fourier space power spectra shown in Figs. [5] and |H] tell 
us the characteristics of the magnetic fields injected into 
the plasma in this model. 

Our results should be viewed as providing initial condi- 
tions for subsequent evolution which will also entail MHD 
and cosmological effects. If the cosmological fluid is tur- 
bulent, say due to the motion of bubble walls, that too 
will play a role. These effects did not enter our study 
because the magnetic fields that are produced due to 
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FIG. 7. Mh(t). Plots are made following the scheme of Fig. [5] 
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FIG. 8. S(k) vs. k. Plots are made following the scheme of 
Fig. El 



sphaleron events are on scales comparable to the inverse 
W-boson mass and far smaller than the scale at which the 
medium can be treated like a fluid. The micrscopic pro- 
duction, however, occurs at a high rate, since it is also 
the rate at which baryons are produced, and the mag- 
netic field due to different sphalcrons will subsequently 
spread, merge and permeate space. 

MHD effects will come into play on length scales that 
are large compared to the thermal scale. As the mag- 
netic field expands to larger scales, but still less than 
the frozen-in scale, MHD effects will change the linear 
expansion to diffusive expansion. On yet larger scales, 
the magnetic field expansion enters the frozen-in regime 
where it can only scale with the expansion of the uni- 
verse. In addition to these considerations, the evolution 
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needs to include the helicity of the magnetic field and any 
turbulence that may accompany the phase transition. It 
is known that helicity can be responsible for an "inverse 
cascade" that transfers power to larger scales. A discus- 
sion of some of these issues in the present context may 
be found in Ref. 0. 

A long term goal of our model for generation of mag- 
netic fields, is to connect the particle physics processes 
during baryogenesis (such as the phase transition) to 
characteristics of the magnetic field. The hope is that 
eventually the observation of a primordial magnetic field 
may say something about particle physics at the baryo- 
genesis scale, the nano-second universe, and perhaps also 
the observed astrophysical magnetic structures. 
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